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Abstract 

Using techniques of complex analysis in an algebraic approach, we solve the 

wave equation for a two-level atom interacting with a monochromatic light 

field exactly. A closed-form expression for the quasi-energies is obtained, 

which shows that the Bloch-Siegert shift is always finite, regardless of whether 

the original or the shifted level spacing is an integral multiple of the driving 

frequency, uj. We also find that the wave functions, though finite when the 

original level spacing is an integral multiple of u, become divergent when 

the intensity-dependent shifted energy spacing is an integral multiple of the 

photon energy. This result provides, for the first time in the literature, an 

ab-initio theoretical explanation for the occurrence of the Freeman resonances 

observed in above-threshold ionization experiments. 
PACS number(s): 32.80.Rm, 42.65.Ky, 12.20Ds, 03.65.Nk 
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I. INTRODUCTION 



The interaction between light and matter is a fundamental problem in physics, whose 
study has led to the birth of quantum theory about a century ago. The two-level atom model 
was originally proposed by Einstein [1] to study the transitions between two energy levels 
of an atom interacting with light (later especially with laser light). However, despite many 
great and significant progresses since then, an exact solution remains elusive even for the 
simplest problem of a two-level atom [2] interacting with a classical or quantum-mechanical 
light field. Pursuing higher accuracy in describing a physical system is always an ultimate 
goal for physicists. 

Exact and analytic expressions also provide new starting points for further developments 
of theories in physics. The exact quasi-energy levels and the wave functions obtained by 
solving the two-level atom model can be used in the calculation of many important physical 
quantities, such as the Rabi flopping frequency and the inversion rate. Mathematically, 
an approach that exactly solves this model may provide a starting point for solving more 
complicated cases, such as a driven N-level atom, as well as an atom in a multi-mode laser 
field. 

Since the pioneering work of Bloch and Siegert (BS) [3], there have been many different 
approximate methods developed to solve a two-level system driven by an external field. 
The rotating wave approximation (RWA) is a widely used method. As is well-known, this 
approximation is good only when the frequency of the light field is "near resonance" but 
also not too close to the resonance. The word "near" means that the frequency of the 
light field is near the original energy spacing of the two-level atom. Many works have been 
devoted to go beyond the RWA method. For example, Shirley [4] applied Floquet's theorem 
and perturbation method to solve the time-dependent Schrodinger equation for a two-level 
system. And Cohen- Tannoudji et al. [5] used perturbation methods to solve a quantum- 
field two-level system. Mittleman et al. [6] have recently derived the wave functions and 
quasi energies for a two-level atom driven by a low-frequency strong laser pulse and applied 
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their result to emission spectra and high harmonic generation (HHG). The low- frequency 
approximation adopted by Mittleman et al. skips all higher resonances and can be thought 
of as a limiting case opposite to the near-resonance approximation. The Continued Fraction 
(CF) method, giving recurrence relations for the Fourier coefficients of the wave function, 
is also a commonly adopted method going beyond the RWA. Swain [7] , Yeh and Stehle [8] , 
Becker [9], and recently Feng [10] applied the CF method to obtain approximate solutions. 
In continued fractions expressing the wave function, the unknown quasi-energy is involved. 
Approximations, used in evaluating the quasi-energy, make the corresponding wave functions 
inaccurate. Due to the infinite order of the algebraic equations satisfied by the quasi energy, 
the exact value of quasi energy in a closed form has never been derived. 

In 1987, R. R. Freeman et al. found experimentally [11] that the above-threshold ion- 
ization (ATI) peaks broke up into many small peaks when the laser pulses were short. 
The appearance of the small peaks were interpreted as multiphoton resonances between 
the ground state, say 5P 3 / 2 for the outermost shell electrons of xenon atoms, and Rydberg 
states with a shifted energy level. In the literature these resonances are now called Freeman 
resonances and their appearance can be phenomenologically explained by a.c. Stark-shifted 
multiphoton resonances. 

Such multiphoton resonances have been observed for years and modelled theoretically 
using Floquet and numerical approaches [12]. Due to interactions with the radiation field, 
the energy level spacing of a two-level atom acquires an intensity dependent shift. In the 
theoretical literature ,the term "near resonance" refers to the condition where the photon 
energy is near the original energy spacing (pre- resonances, in the absence of radiation). How- 
ever, the fine structure of ATI peaks observed by Freeman et al. [11] can be well interpreted 
as Rydberg state resonances occurring only when the shifted energy spacing is equal to an 
integral multiple of the photon energy, the Freeman resonances, in the presence of strong 
radiation fields. It is this intensity-dependent shift in the resonance frequency that calls 
for a fundamental explanation in the theory for a driven atom. Moreover, one is certainly 
tempted to know what happens when the light field is neither very near nor very far away 
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from any resonance. Even in the so-called near-resonance case, calculations and analysis 
with higher accuracy are always desired. All of these require a solution to the two-level 
atom problem that is as exactly as possible. 

In this paper, we attack the two-level atom problem in an algebraic approach. We start 
with a proof of the equivalence between a classical-field description and a corresponding 
quantum-field description for a two-level system in a driving field. Then, we recast the 
classical-field differential equations of motion into an infinite system of linear equations, with 
the energy determinant of infinite rank of the form of a continuant [13], i.e. having non-zero 
elements only on three major diagonals. We then directly evaluate the energy determinant, 
using techniques in complex analysis and the trick that breaks the relevant continuants into 
sub-continuants of a half infinite rank, which are further expressed as infinite series. In this 
way, a closed-form expression for the quasi energies is obtained. 

Our solutions exhibit several interesting features: 1) A simple cosine energy shift formula 
is derived, which naturally exhibits the Floquet quasi-energy feature. 2) It incorporates 
multiphoton effects, in particular all multi-photon resonances if there are any. 3) For the pre- 
resonance cases, when the original energy spacing is an integral multiple of the photon energy 
of the radiation field, both the shifted quasi-energies and the corresponding wave functions 
are finite. Therefore no singularity (or resonance) really appears at the pre-resonances. 4) 
It shows theoretically the existence of Freeman resonances, namely the wave function has a 
singularity when the intensity dependent energy spacing shifts to an integral multiple of the 
external photon energy. 

Comparisons are made between our result and earlier results indeed shows agreement at 
the leading order. The higher order correctness of our results can be guaranteed and checked 
by the mathematical derivation process and also by the comparison with experimental find- 
ings. 
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II. EQUATIONS OF MOTION 



The goal of this paper is to solve the following equations of motion which describes a 
two-level atom driven by a radiation field: 

(-f- -iD cos to x + iAa z + iEI)Y = 0. (1) 

GST 

In this equation, we have chosen c = h = 1, and set the field frequency cu — 1. Here 
2A stands for the energy spacing in units of oo in the absence of the radiation field; we 
have introduced a dimensionless dipole moment D for the interaction strength, with D 2 
proportional to the laser beam intensity. The notation E stands for the quasi-energy, in 
uints of uj, of the two-level atom in the classical radiation field. It can also be directly called 
the energy level, if one treats the radiation as a quantum field. Y = (Y 1 ,Y 2 ) t , and a x , a z 
are Pauli matrices; / is the 2x2 unit matrix. 

This equation is usually derived from the quantum mechanical equations of motion for 
a two-level atom interacting with a classical, single-frequency mode in the dipole approxi- 
mation. Below we will show that it can also be derived from the quantum field approach 
in the large photon number (LPN) limit without any other approximation. For this reason, 
we regard the equations of motion (1) as an exact description of a driven two-level system 
when the driving field is a classical field. 

Let us start from the equations of motion for a two-level atom interacting with a single- 
mode quantun field: 

[Aa z + uN+ \e\ga x {a + <f)]\<j>) = -E\<j>), (2) 

which is equivalent to the one in Cohen- Tannoudji et a/.'s paper [5], if we set 2A = cu , \a = 
= 7^\,E = —E'. 

Now, we introduce a new basis (y\ = J2 n y n ( n \ wrfc h (n\ being Fock states: 

|0)-0(y) = £2/>l</>>- (3) 

n 

Then, in the large-photon-number (LPN) limit, the equations of motion (2) becomes 
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Y\^zV n + ny n cu + \e\AaM~ 1 + = -££z/Vn, (4) 

n n 

where 

0n = (n\4>), A = (5) 



We further rewrite the above equation as 

d 
dy 



[Aa z + uy^- + \e\Aa x (y~ l + y)\<j>{y) = -E(j){y). (6) 



Letting 

2|e|A = Du, y = -e iWT , u = 1, (7) 

we obtain eq. (1). Starting with eq. (1) and going backward through the proof, with 
mapping y n to Fock state \n) and resuming the commutation relation between y and y -1 , 
we can recover eq. (2). To recover ui, one should just take A — > A/c<j and £ — > E'/a;, with 
D remaining dimensionless. 

III. SOLVING THE ALGEBRAIC EQUATIONS OF MOTION 

Write the solutions in the case of A ^ in the form: 

Y = CMr)Y« + C 2 (r)Y( 2 \ (8) 

where Y^ and Y^ 2 ^ are two linearly independent solutions of eq. (1) with A = = E: 

YW(t) = e- i£,sinT 4=(l, -1)*, Y (2) (r) = e iDsinT ^=(l, 1)*. (9) 
v 2 v 2 

Here the superscript t means transposition. These solutions are ortho-normal: YW^Y^' = 

and satisfy 

Y^t^YW = Y^Y^ = 0, 
Y( 1 )t (7 ,Y( 2 ) = (Y^Y^)* = e i2DsiQT 
The resulting differential equations are 



i4-Ci(r) - AC 2 (T)e i2DshlT - Ed(r) = 

GST 

i^-C 2 {r) - Ad(r) e - i2Dsi ^ - EC 2 (r) = 0. (10) 

GST 

Using the following expansions, 

Ci(r) = £ C ls e~^, C 2 (r) = ]T C^e"^, (11) 

Equations (10) can be transformed into a set of linear equations: 

sC ls - A ]T C 2t J t _ a (2£>) - £C ls = 0, 
t 

sC 2s - A C lt J t - a (-2D) - EC 2s = 0, (12) 
t 

where J n (x) are ordinary Bessel functions, s and t are integers running from — oo to oo. 
However, each of the equations involves an infinite sum of terms, so they still look awful. 
To simplify them, we use the Bessel functions to construct a transformation: 

A q = Y,Ci s h- q {-2D). (13) 

s 

The inverse transformation is 

C la = £V*-(2£>)- (14) 

Using the inverse transformation to express eq. (12), and the recurrence relations for the 
Bessel functions nJ n (2D) = D[J n _i(2D) + J n+1 (2D)], we obtain 

AC 2s = -D(A S+1 + A a _x) + (s - E)A S , 

C 2s = j^A s . (15) 
So finally the new variables A s satisfy a set of simple linear equations 

( /_ ( f )2 ~_ S) A2 (^ + i + ^-i) + ^ = 0, (16) 

with s = • • • , —2, —1, 0, 1, 2, • • •. Each equation now involves only three terms. This success 
of simplification is crucial to our subsequent treatments. 



IV. INFINITE DETERMINANT AND QUASI-ENERGIES 



For a nontrivial solution to eqs. (16) to exist, the quasi energy E has to be such that 
the following infinite determinant vanishes: 



a- 2 (3-2( E ) 

7 _!(£) a_x P-!(E) 

7o ( J E?) a A>(£) 

7l (£) «i Pi(E) 

7 2 (£) a 2 



det(E) = 



(17) 



with 



a s = 1, p s (E) = ls+1 (E) = 



D(E-s) 



(18) 



(E — s) 2 — A 2 ' 

A determinant of this type is called tri-diagonal, or it is called a continuant [13]. 

In the previous literature, this kind of infinite determinants were evaluated by using 
various approximations, e.g., the power series expansion in D [10]. Here we will evaluate 
this infinite determinant exactly, using techniques in complex analysis. The key observation 
is that the tri-diagonal infinite determinant in eq. (17) is absolutely convergent, since the 
infinite sum, J2 S Psls+i, is absolutely convergent. (See, for example, the classical treatise 
[15].) 

Therefore, if we regard the energy E as a complex variable, then the infinite determinant 
(17) defines an analytic function on the complex- E plane. Actually it is a meromorphic 
function of E which has two groups of poles at E — ±(A + s), (s — ± 1, ±2, ■ • ■). 

When 2 A = n for integer n > 1, we call the case as pre-resonance. In the non-pre- 
resonant case where 2 A ^ n, the poles of det(E) are all simple poles. Since (5 S {E) and 
1s(E) are periodic functions of E with period unity. The poles in each group, E = ±(A + s) 
respectively, are equally spaced and have the same residue. The residues, 
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R+= lim (E — A±n) det(E), R- = lim (E + A ± n) det(E), (19) 

are independent of n. Furthermore, it can be directly verified that the residues R + and _R_ 
differ from each other only by a sign: 

R + (D,A) = -R_(D,A). (20) 

Each group of poles with the same residue on the real E axis of the complex-E plane 
suggests a cotangent function. We also observe that det(E = ioo) = 1. Thus the exact value 
of the determinant (17) can only be the following function: 

OO 1 1 

det(E)= 1 + R+ V (- - A - A ) 

v ; + ^ K E-A + n E + A + n 



n=—oo 



= 1 + tiR + cot[7i(E - A)] - nR + cot[7i(E + A)] 
sin(2vrA) 

+ + sin[7r(E - A)] sin[7r(£ + A)] ' ( } 

By the uniqueness theorem in complex analysis, the right side of eq. (21) and that of eq. 
(17) agree with each other on the whole E-plane. 

When 2A is a positive integer, the two groups of simple poles are merged to become 
double poles. When A = 0, the poles at E — integer still remain as a first-order or zeroth- 
order ones. (The details will be discussed in Appendix.) 

Thus we are able to put the characteristic equation, det(E') = 0, into an exact and closed 
form: 

sin[7r(£ - A)] sin[7r(£ + A)] = -ttR + sin(27rA) . (22) 

Solving the above equation, with inclusion of the pre-resonance case, we obtain a cosine 
function of quasi-energy: 

cos(2tt£) = cos(2ttA) + 2nR + sin(27rA), (2A ^ n) 

cos(2tt£) = cos(2ttA) + (-l) n 27r 2 r n , (2A = n), (23) 

where r n are residues of R + as a function of 2A, defined by the following limiting processes: 



lim #+sin(27rA) = (-1)" lim R + n(2A - n) = (-l) n nr n . 

2A^n 2A-+n 



(24) 



We note that r = for n = 0. (The detailed proof is presented in the Appendix.) 

This expression (23) for energy shift has a unique advantage: the Floquet condition is 
automatically satisfied in view of the cosine function. 

The next main job is to evaluate the factor R + . In the context below, for writing 
convenience, we use finite determinant notations to express the infinite determinants. We 
find 



Jiin (£ - A) det(E) 



D(A+2) 
2(2A+2) 



D(A+3) 
3(2A+3) 



D(A+1) 
(2A+1) 



D(A+2) 
2(2A+2) 



D 
2 



£>(A+1) 
(2A+1) 



J>(A-1) 
-1-(2A-1) 



D 
2 

1 

D(A-2) 
-2(2A-2) 



P(A-l) 
-1-(2A-1) 

1 

D(A-3) 



D(A-2) 
-2(2A-2) 

1 



(25) 



-3(2A-3) 

The following lemma is useful to determine the residues of R + : 

LEMMA: For any integer n, the residues of det(-E) with 2 A = n at the second order 
poles E = l + n/2, where / is an arbitrary integer, are equal to the residue of R + (2A) at the 
first order pole 2 A = n; i.e., 



E— >(Z+n/2) 

where r„ are residue of the function R + (2A) at the pole 2 A = n. 
Proof. From eqs. (21) and (24) 

(-l) n 7r 2 r n 



(26) 



det(£)| 2A= 



1 + 



sm[n(E - A)]sm[Ti(E + A)] l2A - n " ' sin 2 [7r(E - I - n/2)} ' 

we can see that when £'—>(/' + n/2), the function n 2 / sin 2 [7r(E — I — n/2)] behaves like a 
second order pole [E — (I + n/2)] -2 . 



1 + 



nr. 



n 



(27) 
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In the n = case, the above proof still holds, since -R+|2A=o = (see appendix) and 
r = 0. QED 

From this lemma we learn that the function of R + (2A) can only have first-order poles 
at non-zero integers. Thus, R + can be expressed as 

4A 



R + = Rr+ e 



E 



(28) 



, _ X) 2A-n £i4A'-n* 

where (a) r n = r_ n , (b) r = 0, and (c) R+ = have been used. The proofs for (b) and (c) 
are given in Appendix. The proof for (a) is the following: Using eq. (26) to express r_ n , we 
verify that the values of the two factors on the left hand side of this equation do not change 
with changing n — > —n. That the second factor does not change can be seen from eqs. (21) 
and (20) while the first one can be seen with substituting I in I — n/2 by I' + n since both I 
and I' can be arbitrary integers. 

Further evaluation of R + is based on the evaluation of the residues r n . In the following, 
we list few low-photon-number residues as examples: 



1 ¥t 



3-4 



2-3 



D 



5 
_2 
2-3 



D 



3 
2. 
1-2 



D 



3 
2. 
1-2 



D 
2 



D 
2 



ri 



D 
2 



D 

2 



(29) 



D 



3 

2-1 



¥t 1 -¥t 



D 



5 
2. 
3-2 



2-1 



4-3 



5 

2. 

3-2 



11 



r 2 



3D 
2-4 



i 4D 
3-5 



2D 
1-3 



3D 
2-4 



D 
2 



2D 
1-3 






D 
2 



D 
2 







2D 
31 



D 
2 

1 

3D 
4-2 



_4D I 
5-3 



(30) 



2D 
31 

1 



3D 
4-2 



?~3 



and 



r 4 



1 ^ 



D 2 

5-2 



D- 

^2 I " 2 



I D4 

6-2 



D 
2 



Dj 
4-1 



^ 

2 u 



2-(-l) 



D3 I D3 
5-1 5-1 



Dl 



3-(-l) 



D 
2 



1 

l.(-2) 



D 
2 

1 





Dt 



2-(-l) 
1 



Dl 
3-(-l) 

1 

J>(-1) 
l-(-3) 



D 
2 





1 

D 
2 



l-(-2) 


(-l)(-4) 



l-(-3) 


g(-3) 
(-l)-(-5) 



D 
2 

1 

(-2)(-5) 



D 
2 

1 

DM 



(-l)(-4) 

1 



(31) 



g(-3) 
(-l)-(-5) 

1 



(-2)-(-6) 

For general n, r n can be expressed as an determinant of infinite rank 



(32) 



12 



£>(f +2) 
2(n+2) 



g(§+3) 
3(n+3) 

1 

l(n+l) 



D(f+2) 
2(n+2) 

1 

_D 

2 



g(f+l) 
l(n+l) 



-l(n-l) 



2 



g(§-2) 
-2(n-2) 



-l(n-l) 
1 

-3(n-3) 



-2(n-2) 
1 



1 

g(f+2) 
2(n+2) 



g(§+3) 
3(n+3) 

1 

l(n+l) 



D(f+2) 
2(n+2) 

1 

V 
2 



l(n+l) 


-l(n-l) 



2 
1 

g(j-2) 
-2(n-2) 



-l(n-l) 
1 

-3(n-3) 



J>(?-2) 
-2(n-2) 

1 



(33) 



In this notation, the lower part in the expression is the continuation of the upper part in 
the direction of the main diagonal. 

We have been able to work out the first a few terms for the residues in the dipole 
expansion. Here we only cite the results: 



ri 



4 1 48 64 ) 



(n^l), 



(34) 



The numerical coefficients given here are all exact. 
Thus, the expression for R + up to the D 4 term is 
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4 48 64 



2 1 \ 

+ 777 +••• 



4A ~ 



D 4 



4 \n 



n 



+ .. 



4A 



4A 2 - n 2 



(35) 



The explicit form of R + can be used to evaluate the energy shift. In the pre- resonance case 



cos(2tt£) 



7T 



7T 



24 32 



cos(2tt£) = (-1)"- (-1) ; 



2 W 



n 



(2A = 1) 
(2A = n^ 1). 



(36) 



There are many exact ways to express E as an exact function of A. Following are the 
suggested ones: 



1 



E = — cos' 1 Jcos 2 (ttA) + ttR + sin(27rA), (2A ^ n) 
ix v 

E = -cos-^TrV-rafe+i) (2A = 2fc + l), 

7T 



(37) 



and 



£ = i sin -1 ^/sin 2 (7rA) - 7ri? + sin(27rA), (2A ^ n) 

(2A = 2ifc). 



7T 
1 



E = - sin 1 (7r v /= r^ 

7T 



(38) 



The right side of the equations can all be added with an integer, due to the Floquet feature 
for quasi-energy. 

The Bloch-Siegert shift, defined as E BS = 2E — 2A, has now the exact expressions 



E BS = - cos" 1 t/cos 2 (7rA) + uR+ sin(27rA) - 2A, (2A ^ n) 
E BS = -cos-\n^F^)-2A (2A = 2/c + l), 

71 



(39) 



and 



E BS = - sin" 1 ^/sin 2 (7rA) - 7ri?+ sin(27rA) - 2 A, (2A ^ n) 



7T 

2 . , 

E BS = - sin-^Trv^r^) - 2A (2 A = 2fc). 

7T 



(40) 
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V. THE WAVE FUNCTIONS 



The equation (16) can be further written as [14] 
for s > 0, 



A s -1 



for s < 0, 



A s ^ (E - s) 2 - A 2 A s+1 ' 

- s) A s 
A„ -1 



A a+1 (E-s) 2 -A 2 ' 
D(E-s) + A s 



A_i £ 2 -A 2 A 1 

for s = ' ^ + ^r + ^ = a 

By iterating the first relation, we express Ai/A as a continued fraction: 

A 1 

A (E-l) 2 -A 2 1 



D(E-l) (E-2) 2 -A< 



D(E-2) 

Similarly by iterating the second relation, we have 
A_i 1 



(E + l) 2 — A 2 1 

~ D(¥TT) (E + 2) 2 - A 2 



£>(£ + 2) 

Putting back all the transformations made, we obtain 

YT = EE ^J 9 _ s (2L>)e-^Y( 1 ) + £ ^^e— Y( 2 ) 
= E ^e-* aT e i2J?sinT Y( 1 ) + 2 7^^ e "" TY(2) 

= J2 ^ s T e^ ST e i2Dsinr e- lDsinr ^(l, -1)* + J2 ^TT^e- iST e lDsinT ^(l, 
s v2 s s±E V2 

= £ ^-'"-^[(1, -1)' + ^(1, l)<] e <°*- 
where the superscript =f denotes the solutions corresponds to ±E. 
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VI. FREEMAN RESONANCES 



In our solutions, we identify a Freeman resonance when the new energy spacing 2E, which 
is field-intensity dependent, is an integral multiple of the field photon energy The derived 
wave functions apparently have singularities only at 2E = 2shw (s = 0, ±1, ±2, ±3, ...), and 
not at the pre- resonance case. Here we do not exclude the 2E = case, since 2E = means 
2E can be any integer. We also see from the wave function that in the pre-resonance case, 
the wave functions are finite, as well as the quasi energies given by eqs. (37) and (38). 

All this means that the pre-resonances are not true resonances; only Freeman resonances 
are true resonances. At first glance, one may think that the Freeman resonances occur only 
when the resonating photon number is an even number. Since the quasi-energy spacing 
2E can be added with an arbitrary integral multiple of the photon energy, we can replace 
2E' = 2E + 1 in above equations related to the wave functions. Thus, we immediately 
find the resonances with odd photon numbers. This analysis also indicates that we need 
four or more different quasi-energy levels as basic ones even in the non-resonance case, since 
2E' = 2E + 1 may give different wave functions. On the other hand, we do not need more 
than four as basic ones, since 2E' = 2E + 2 will not give a new wave function and it is 
included in the iteration process in the continued fractions. Thus, we conclude that we need 
four and only four different quasi-energies as basic ones to produce the wave functions. 

At the Freeman resonances, the intensity of the field and the original energy spacing 2A 
satisfy the following equation from eq. (23), 



cos(2ttA) + 2nR + sin(27rA) = (-l) n . 



(47) 



When n = odd, the equation reduces to 



2irR+ = -cot 2 (yrA). 



(48) 



When n = even, the equation reduces to 



2%R + = tan 2 (7rA). 



(49) 
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In this case, as we pointed before, the wave functions have an infinite amplitude. 

Equations (48) and (49) are transcendental equations. For a given field intensity and 
a Freeman resonance (E = rtw), one can solve the transcendental equation to identify the 
resonating atomic level, which has the original spacing 2A from the ground state. This 
equation has discrete solutions for 2A, which change when the field intensity changes. This 
theoretical feature does agree with experimental findings. Experimentalists call an electron 
energy peak in an ATI spectrum a Freeman resonance when the energy spacing between 
the ground state and the energy peak is an integral multiple of the laser photon energy, 
and intepret the energy peak as a formation of photoelectrons from a shifted Rydberg state. 
ATI spectra do show that the Freeman resonances have a discrete feature for a fixed laser 
intensity. When the laser intensity changes, different sets of Rydberg levels come into play 
consecutively as Freeman resonances appearing in the ATI spectrum. 

VII. DISCUSSIONS 

The basic requirement to a correct solution of an interacting system is when the interac- 
tion vanishes, the solution reduces to the one of the corresponding non-interacting system. 
For the problem in hand, the interaction is imposed through the dipole moment, D. Thus 
the leading term of quasi-energy in the expansion in powers of D 2 u 2 should satisfy the basic 
requirement, also should signify the physical meaning of field intensity. Another important 
quantity is (2 A — uj) 2 which signifies the detuning of the field frequency from the transition 
frequency of the original two-level atom and competes with D 2 oj 2 in the near pre-resonance 
processes. To see the interesting competition in different limiting processes, we consider the 
following two cases respectively. 

A. Energy Shift at Exact n-photon Pre-Resonances 

Practically, if we can tune up a laser in a way that an integral multiple of the laser 
frequency, nuj, matches the energy spacing, 2 A, of a two level atom i.e., (2A — nuo) 2 <C D 2 uj 2 , 
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we can set 2 A = nuo and use eqs. (36), (37), (38), (39), and (40) to obtain the energy shift 
in the small-dipole limit. 

We have the following three sub-cases: 

I. The Single-Photon Case 

In the single-photon pre-resonance case, the quasi-energy E has a simple form by just 
keeping the leading term of D, with u recovered in the expression explicitly, 



£ = r os l±- )■ (50) 

Using cos _1 (x + d) = cos -1 x — (1 — :r 2 )~ 1//2 <i, where x = and d = ±%D/2, we find 

E = ^(l + D), E bs = ujD. (51) 

where the negative sign in eq. (50) is selected according to the limiting process from the 
near pre-resonance case in the next subsection. 

This is a simple interesting results that the shift is proportional to D. 

2. The Odd-Number Photon Case (n ^ 1) 

E=-cos- 1 [± . 52 

Using cos _1 (:r + d) — cos -1 x — (1 — x 2 )~ 1 / 2 d, where x = and d = ±nnD 2 /2(n 2 — 1), we 
find 

with the same sign selection rule as above. 

3. The Even-Number Photon Case 



£ = - sin 1 ± . (54) 

7T I 2 n 2 - 1 ' 



Using sin l (x + d) = sin ^ + (1 - x 2 ) l/2 d, where x = and d = ±nnD 2 /2(n 2 - 1), we 
find 

nu . D 2 , ^ ncj£> 2 . . 

B = lT (1 + rf £ - = ^T- (55) 



where the positive sign is selected in eq. (54) according to the same limiting process as 
dsicussed above. 

Now we see that the BS shift is proportional to D in the single-photon pre-resonance 
case, while it is proportional to D 2 in the multiphoton pre-resonance case. 

B. Energy Shift in the Weak-Field and Near Pre-Resonance Case 

In the previous case, we let (2A— nuj) 2 be infinitesimal first. In the rpesent case, switching 
the limiting procedures we let D 2 u 2 be infinitesimal first, i.e., D 2 uo 2 <C (2A — nuo) 2 . 
To treat this case, we need to expand the term 2ir sin(27rA)_R + in eq. (23) as 

1 1 , 



27rsin(27rA) ^(-—— + 



r.; 



f2A-j ' 2A • ./"•' 



1 ^ 1 



(-ir2 7 r 2 {r ri + (2A-n)( £ — - + £ 



3 



r 



3 



2 i oo i 

-<2A - „) V . - (2A - „)>[£ — 2 + g (56, 

We consider the following sub-cases: 
1. The Single-Photon Case 

In the small-dipole case we only keep r 1; because only r\ has D 2 as its leading order, 
while all other r n have D A as their leading order. Thus, from eqs. (23) and (28), the cosine 
energy relations read 

cos(2tt-) = -1 + ^(— - l) 2 - 2n 2 n - n 2 n (— - 1), (57) 

where ri —D 2 /4 in the leading order. By expanding cos(2ttE/uj) ps — 1+tt 2 [(2E/uj) — l] 2 /2, 
we obtain an approximated quadratic equation 

D 2 

(2E - uj) 2 = (2A - uj) 2 + D 2 uj 2 + — cj(2A - uj), (58) 

2 

which has solution 



D 2 

Ebs = {2E - 2A) = -(2A - u) ± W(2A - uj) 2 + D 2 uo 2 + — cj(2A - (59) 
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This result agrees with the one from RWA but with an extra correction term. 
For a small dipole D 2 <C 1, choosing the positive sign, the BS shift reduces to 

The extra term signifies the one-photon emission effect beyond the RWA. Here we see that 
when D 2 uj 2 <C (2 A — uj) 2 , i.e., the interaction term is smaller than the detuning, the BS 
shift is proportional to D 2 . 

The sign selection here is forced by the requirement that when the interaction Duo — > 
0, the value of E BS has to vanish. With the positive sign selected in eq. (59), letting 
(2 A — u) — > directly, we get the same result as eq. (51), whose sign is thus determined. 

2. The Multiphoton Case 

We have the following expansion: 

cos(2tt-) = (-1)»{1 - ^(— - n) 2 } + (-lT2n 2 (^— + -^-)(— - n)n. (61) 
to 2 uj n — 1 n + 1 uo 

By expanding cos(2nE/uo) xs (— 1)™{1 — tc 2 [(2E/u) —n] 2 /2}, we obtain an approximated 
quadratic equation: 

2n 

(2E - nu) 2 = (2A - nuj) 2 + D 2 u— (2A - nu), (62) 

n 2 — 1 

which has teh solution 

E BS = (2E - 2A) = D 2 lu^—, (63) 

where the sign is determined by the same method as before, which also justifies the sign 
selection in eqs. (53) and (55). An amazing thing here is that in the two different limiting 
processes and with different formulas, Ebs behaves the same way, both proportional to D 2 
and with the same proportionality constant! 

VIII. CONCLUSIONS 

In the following we summarize the features of the exact solution obtained: 
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1) It explicitly exhibits the Floquet quasi energy behavior, namely the quasi-energy 
spacing is determined only up to its cosine. 

2) It incorporates all multi-photon effects. 

3) When the original energy spacing matches an integer number of photon energy, the 
pre-resonance case, the interacting system does not really resonate. Both quasi-energies and 
wave functions are finite. 

4) We have shown theoretically that when the shifted energy spacing matches an integer 
number of photon energy, the interacting system resonates. This theoretical feature explains 
the Freeman resonances observed in ATI experiments. We have also proven that the Freeman 
resonances have a discrete feature, and when the field intensity changes different sets of 
unoccupied, excited atomic levels come into play as the resonances. 

5) The Bloch-Siegert shift vanishes when the interaction D vanishes. The way it vanishes 
obeys some simple rules: Near the single-photon pre-resonance, whether the BS shift vanishes 
as D or D 2 depends on whether the detuning first goes to zero or the interaction strength 
first goes to zero. Near the multiphoton pre-resonances, the BS shift always vanishes as D 2 , 
independent of the limiting process. These rules can be subject to experimental tests. 

Acknowledgments DSG was supported in part by a summer grant of Physics Depart- 
ment of Ohio State University. YSW was supported in part by the NSF through grant 
PHY-9970701 and PHY-0457018. 

IX. APPENDIX 

We will show that R + as the residue of det(E') vanishes when 2 A = 0, and that 
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R™ = lim R+ 
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In the case of 2 A = 0, 
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D 
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We notice 



(64) 



(65) 



hm(E-n) det(E)\ A=0 = R + (A = 0), 

E^n 



(66) 



which does not depend on n; thus, 



det(E)\ A=0 = 1 + 7ri? + (A = 0) cot(nE). 



(67) 



The characteristic equation is obtained by setting det(E)\ A=0 — 0, i- e -> 



sin(7r£) + nR + (A = 0) cos(tt£) = 0, 



(68) 



which has solutions 



E = tan^ 1 (7ri?+(A = 0)), 

7T 



(69) 



where R+(A = 0) is 
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1 # 
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D 
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The last step ( = ) needs a proof. 
Proof: 

We use a recurrence relation of Bessel functions 



x 



D 
' 1 



¥ 1 



D 
' 3 



7^— (Jn-1 + Jn+l) + Jn — 0, 

In 

K{J. X + Ji) = 0, 



= 0. 



" 2 



(70) 



(71) 



where X is an arbitrary constant. 

Let n run from — oo to oo. The algebraic equation set for J n has non-zero solutions. So 
the coefficient determinant must be zero, i.e., 



det(x) 



1 if 








¥ 1 


— X 

4 






— X 
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— X 

2 






K 


if • • 


= 



1 2 . 

2 

4 1 4 



(72) 



By setting x = 
get eq. (70). 
QED 



-D and K = D/2we get eq. (64). By setting x = -2D and K = D/2we 
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